Effect of doping on irreversible Campbell length in MgB 2 single crystals 
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We have measured small-amplitude rf penetration depth \{H, T) in pure and C, Li and (Li+C) 
doped single crystals of MgB2 . The effect of doping on the critical temperature T c and on the upper 
critical field H C 2 was found to be in good agreement with previous results. We report the presence 
of clear signatures of irreversibility in \{H,T), associated with the peak effect. Carbon doping 
enhances the signature and shifts its position on the H-T phase diagram to higher temperatures. 
Contrary to this, Li substitution has the effect of suppressing the peak effect to lower temperatures. 
Analysis of both zero field (ZFC) and field cooled (FC) measurements suggests that the hysteresis 
associated with the peak effect is due to macroscopic screening supercurrents, j, generated during 
the magnetic field ramp. 

PACS numbers: 74.25.-q,74.25.Qt 



o 
o 

I 

u 

Oh 

=5 

CZ2 
+-* 



I 

C 

o 
o 



> 

00 

o 
en 

i> 

o 

oo 
O 

> 

X 

S-H 



I. INTRODUCTION 

MgB2 is a two-band (two quasi 2D, er-bands and two 
3D, 7r-bands), two-gap (c-gap, A CT w 78 K and 7r-gap, 
A,,- rs 29 K) superconductor with high critical temper- 
ature T c « 39 K for an s-wave superconductor 1 . These 
characteristics have made MgB 2 very attractive for ad- 
vancing theoretical understanding of superconductivity 
in two-gap systems. Theory has shown that the two gap 
nature has a significant effect on the H C 2(T) phase dia- 
gram, which is governed by the ratio between the intra- 
band diffusivities D a /D^. The high transition temper- 
ature is advantageous for practical applications, however 
the potential is limited by the relatively low values for 
the upper critical fields, particularly for the c-axis orien- 
tation, H^ = 3-6 T. 

Atomic substitution is a method to affect different 
physical quantities like carrier concentration, density of 
states, lattice parameters, interband and/or intraband 
scattering rates in a controlled way. Examining how 
these quantities affect the superconducting properties is 
important both for theoretical understanding and for im- 
proving the potential for practical use of MgB2- More- 
over, doping may influence significantly the vortex lattice 
st ruc ture an d dynamics. Successful substitutions with 
APU, MiPSJand LP for Mg, and C01 for B, as well as 
co-doping d 10 1 11 1 have been reported. 

Hall probe measurements of local ac-magnetic suscep- 
tibility revealed the existence of the peak effect (PE), 
i.e. a maximum of the diamagnetic screening in a cer- 
tain region of the H-T phase diagram, in single crystals 
of both pure and C-doped and Al-doped MgB j 12 ' 13 1 1 ^. 
This effect had been previously observed in both high-T c 
and low-T c superconductors, but theoretically its origin 
is not yet fully understood and several possible expla- 
nations have been put forth 15 . Larkin and Ovchinnikov 
proposed that the effect is caused by a softening of the 
elastic moduli of the flux-line lattice (FLL) in the vicin- 



ity of H c p2l. Alternatively, it was suggested that as the 
elastic moduli softens, the pinning forces overcome the 
elastic ones, marking a transition from weak to strong 
pinninglSl. 

A completely different interpretation was generated by 
the strong experimental e viden ce for the peak effect in 
disordered systems, like NrJ 19 * 20 l Earlier theoretical work 
on vortex matter in superconductors with disorder in- 
duced pinning had claimed that instead of the Abrikosov 
lattice, i.e. long range ordered vortex state, a quasi- 
ordered, Bragg glass phase stabilizes^. Therefore the 
peak effect feature was interpreted as a transition from 
Bragg-glass to a disordered vortex phase. The hysteresis 
of ac-susceptibility which develops at the PE, particu- 
larly the low diamagnetic screening state which results 
after zero field cooling (ZFC) the sample and applying 
magnetic field, was regard ed as evidence for the Bragg 
glass vortex phase in MgBjpnSl. A recent experimental 
study^ using the same technique involved in the present 
work found similar hysteresic behavior of the rf pene- 
tration depth in a series of high and low temperature 
superconductors. However, Ref. [TB] suggests that ramp- 
ing the magnetic field after ZFC gives rise to macroscopic 
screening supercurrents, j, which shift the vortices into a 
state of inhomogeneous distribution, in accordance with 
the critical state (Bean) model. The resulting state is 
a displaced vortex lattice that is thought to disappear 
when field cooling (FC) the sample due to a relaxation 
of screening currents, therefore giving rise to hysteresis. 

The present work studies the effect of C, Li and (Li+C) 
doping on the vortex dynamics in single crystals of MgB2. 
The results of our study are two-fold: (1) We show that 
the same anomaly of rf susceptibility described in the 
previous paragraph is present in both pure and doped 
MgB 2 and it coincides on the H-T phase diagram with 
the measured peak effect from local ac-ma gnetic suscep- 
tibility from Hall probe experiment s ' 12 * 13 ! (2) We show 
that the C doping at B site enhances the maximum of 



the rf-diamagnetic screening by increasing the number of 
pinning centers. Contrary to this, when Li is substituted 
for Mg, the peak effect is suppressed to lower tempera- 
tures, suggesting a slight reduction of pinning. Further, 
studies with both FC and ZFC are shown to be in good 
agreement with the proposed scenario of relaxing super- 
currents after field ramp. The model from Ref. [16] is 
discussed in the context of collective pinning theory as a 
possible explanation for both the shape of the anomaly 
and the history dependence associated with the peak ef- 
fect. 



II. EXPERIMENTAL DETAILS 

The samples used for the present study were single 
crystals of pristine MgB 2 , 12% Li substituted at Mg 
site (Mgo.8sLio.12B2), 4.8% and 7.5% C substituted at 
B site (MgB1.904C0.096 and MgBi. 85 C .i5 respectively), 
and 12% Li and 6% C co-doped (Mgo.8sLio.12B1.8sCo. 12). 
The crystals were grown using a high pressure cubic anvil 
technique, described in Refs. |HI22j . The current sam- 
ples come from the same batches as those studied in 
Ref. [11] where detailed information about the crystal 
structure and stoichiometry are provided. 

The rf susceptibility was measured by placing the sam- 
ple inside the coil of a self-resonant LC circuit, powered 
by a tunnel diode. The resonant frequency of the Tun- 
nel Diode Resonator (TDR), 27r/ = 1/VLC is about 
14 MHz. The experimental data were obtained using a 
3 He refrigerator with a 90 kOe superconducting magnet. 
For the present work, both the dc and the rf magnetic 
fields were applied along the crystallographic c-axis, so 
that the penetration depth was measured in the a6-planc 
(A„ 6 ). 

When a superconducting material is placed inside the 
inductor it produces a change in resonant frequency pro- 
portional to the susceptibility, and in turn, to the pen- 
etration depth of the sample, A/ ex A A. Given that 
the sample shapes were almost rectangular, following the 
procedure described in Ref. [23], the rf susceptibility is 
obtained from: 
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where N is the demagnetization factor and R is the ef- 
fective dimension of the sample as described by Ref. [23] . 
At low temperatures and for applied magnetic fields 
H^> H c i, the penetration depth, and thus the dynamic 
susceptibility, is A 2 = A| + A 2 , where \l is the London 
penetration depth and A„ is the vortex contribution. In 
the pinning regime, the main contribution to the total 
penetration depth is given by the vortex motion, which in 
turn reduces to the Campbell pinning penetration depth 
Ac, resulting from the oscillating Lorentz force exerted 
by the rf screening currents on the vortices: 



X 2 
A C 






(2) 



where 4>$ is the flux quantum and k p is the Labusch pa- 
rameter, which measures the curvature of the pinning 
potential, V(r), n p = \d 2 V(r)/dr 2 \. 

It should be mentioned that, besides technical details, 
there are two major differences between the TDR and 
other ac techniques, such as the mutual inductance or 
based on Hall sensor. One is the frequency range; the 
mutual inductance technique uses low frequencies, up to 
kHz, whereas the TDR is used in our case at MHz val- 
ues. However, this difference should not be important as 
long as both frequencies are significantly lower than the 
pinning frequency, most commonly in the GHz range. A 
more important difference is the magnitude of the exci- 
tation ac magnetic field. As it can be seen from Eq.[2] by 
shaking the vortices with a probe ac magnetic field, i.e. 
ac screening current, the curvature of the pinning poten- 
tial at the vortex location can be determined. The ex- 
cursion that the vortices make while oscillating depends 
on the magnitude of the ac magnetic field. In our exper- 
imental set-up, the applied rf field was estimated to be 5 
mOe < b r f < 10 mOe. For comparison, the ac mag netic 
fields used in two studies of the peak effect in MgBji 2 ^^ 
are larger than 1 Oe, which is at least 100 times larger 
than in the current work. The implication is that the am- 
plitude of the vortex oscillation around the equilibrium 
position are at least 100 times larger in Refs. [12I13J . thus 
the curvature of the pinning potential is not probed lo- 
cally, but rather is averaged over the length of the vortex 
excursion. 



III. RESULTS AND DISCUSSION 

Figure [T] shows the dynamic magnetic susceptibility for 
12% Li substituted (a) and undoped (b) MgB2. It can 
be seen on both figures that for applied magnetic fields 
HoH > 12 kOe, there is an anomaly, i.e. a small kink in 
x{T). For the pristine MgB2, the position of the kink, 
(Hp,T p ), on the H-T phase diagram is in good agree- 
ment with the observed peak effect from Ref. [12] (see 
inset of Fig. [3]). When 12% of Li is substituted for Mg, 
the critical temperature decreases to T c «35K, as previ- 
ously reported in Ref. |llj . However, it is evident from 
Fig. Ilia) that the size of the anomaly is not significantly 
affected by doping. The normalized H-T phase diagram 
in Fig. [4] shows that the PE is suppressed to lower tem- 
peratures in the Li-doped crystal. A possible explanation 
is that doping with Li may actually reduce disorder by 
filling Mg vacancies possibly present in pure MgB 2 crys- 
tals, which are believed to be magnesium deficient. As a 
consequence, the pinning of vortices in the Mg layers is 
reduced, and lower energy thermal fluctuations overcome 
the pinning forces, smearing the PE. This is supported by 
crystal structure results, where a narrowing of some x-ray 
reflection peaks was found in lithium doped samples^. 
Also, the upper critical field H c2 (0) is slightly lower in Li 
doped crystals indicating less disorder. Further de Haas 
-van Alphen measurements to determine the mean free 
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FIG. 1: (a) The rf susceptibility x(r) of Mgo.g8Lio.12B2 in ap- 
plied magnetic fields from 2 kOe to 54 kOe in steps of 2 kOe. 
The position of the peak effect is marked as T p . (b) The rf 
susceptibility x(T) of the undoped MgB2 in magnetic field 
from 14 kOe to 52 kOe in 2 kOe steps. Inset: Zero field x( T ) 
showing the entire transition for the studied samples: MgB2 
(circles, T c «38.5K), MgB1.904C0.096 (crosses, T c «35.5K), 
Mgo.88Lio.12B2 (diamonds, T c «35K), MgBi. 85 Co.i5 (tri- 
angles, T c ^32. 8K), and Mgo.88Lio.12B1.8sCo. 12 (squares, 
T c ^32K). 



path would help with clarifying this hypothesis. 

Contrary to the Li doping, when C is substituted for B, 
the effect on the vortex dynamics is significantly stronger, 
as shown in Fig.pTa). A sharper, more pronounced max- 
imum of the diamagnetic screening develops, similar to 
previous observations from local ac susceptibility at lower 
frequency^. We now compare two different C concen- 
trations: 4.8% (inset of Fig. ga)) and 7.5% (Fig. ga)). 
First, we note that at approximately the same field where 
the PE appears, there are two transitions at T c visible 
in both samples: a lower one associated with the onset 
of superconductivity T C 2, and an upper transition tem- 
perature which we define as T C 3, quite possibly due to 
surface superconductivity. The evolution of T c2 with ap- 
plied magnetic field agrees well with that from resistivity 



FIG. 2: (a) \(T) of MgBr.g5C0.15 at magnetic fields from 10 
kOe to 86 kOe in 2 kOe steps. Inset: x(T) of MgB1.904C0.096 
at magnetic fields from 10 kOe to 58 kOe in 2 kOe steps, (b) 
X(T) for the co-doped sample, Mgo.8sLio.12B1.88Co. 12 in fields 
between 2 and 78 kOe, in steps of 2 kOe, showing very similar 
behavior like the C doping. The traces at 46 kOe (and at 40 
kOe for the inset) have been highlighted for more clarity of 
the peak effect. 



measurements^. Figure [I] suggests that the position of 
the peak effect shifts to higher magnetic fields as the C 
content increases and it appears at significantly larger 
fields than in the undoped or Li-doped MgB 2 . 

For clarity, we have plotted separately, in Fig. [3] the 
phase diagram for the two C-doped samples. It is impor- 
tant to notice that the three main curves, H p (T), H C 2(T) 
and H c3 (T) do not appear to meet at a tricritical point 
in our data both for 4.8% and 7.5% C doping. This is 
in contrast with the results obtained in Nb single crys- 
tal^ 20 , despite the strong similarity between the magnetic 
ac-susceptibility of MgB 2 and Nb. Our results support 
the interpretation that the peak effect is suppressed by 
thermal fluctuations®! and not necessarily at a tricritical 
point where a possible Bragg glass disordering transi- 
tion occurs. Moreover, closer examination of our data 
in Fig. [31 reveals that the PE becomes more pronounced 
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FIG. 3: The phase diagram for the C-doped samples: 
MgB1.g5C0.15 (main graph) and MgB1.904C0.096 (inset). The 
three main curves representing H C 2(T), H C 3(T) and H P (T) are 
shown in both graphs. 



FIG. 4: The reduced phase diagram for the studied samples, 
showing the values for H C 2(T) and the peak effect H P (T), 
normalized to H C 2(0) and to T c . The symbols are explained 
in the legend. Inset: The absolute values of H C 2(T) and H P (T) 
for the undoped MgB2. The values for the peak effect H P (T) 
from Ref. [12: are displayed for comparison. 



and survives to higher temperatures in the sample with 
the higher concentration of C. We attribute this to the 
increase in pinning strength as carbon content increases. 

The susceptibility of the co-doped sample, 
Mgo.88Lio.12B1.88Co. 12, is consistent with the cumu- 
lative effect of Li and C (Fig. pjb)). The peak effect 
is marked by a pronounced dip of x(T) and there are 
two transitions visible at T c , similar to what is seen in 
the carbon doped samples (Fig. [2Fb)). Furthermore, 
the presence of Li shifts the position of the peak effect 
to lower temperatures on the reduced phase diagram 
in Fig. El The co-doping strongly reduces T c (inset of 
Fig. flTby) . A detailed analysis and possible explanation 
for the evolution of T c with doping is given in Ref. p} . 

An important characteristic of the peak effect is the 
presence of hysteresis in dc magnetization and ac suscep- 
tibility. In order to verify and understand the nature of 
the hysteresis, we performed measurements of rf suscep- 
tibility both by zero field cooling- field warming (ZFC- 
FW) and field cooling-field warming (FC-FW) the sam- 
ples, following a scenario schematically illustrated in the 
upper inset of Fig. [5] The results for both C-doped crys- 
tals are shown in Fig. [5] After ZFC the sample to base 
temperature T <C T c , the magnetic field was ramped to 
a certain value resulting in state 1. Then, as the tem- 
perature is increased on curve 1 to 2, the susceptibility 
shows a diamagnetic peak at T p , before the transition 
to the normal state at T C 2 or T C 3. Upon FC the sam- 
ple from 2 to 3, the strong peak is no longer present, 
but at the same temperature T p , the susceptibility devel- 
ops strong hysteresis, resulting in a stronger diamagnetic 
screening at base temperature. By further increasing the 



temperature from 3 to 4, 
most the same path as 2- 



the susceptibility follows al- 
>3, showing only a weak kink 



at T p . Cooling back in field, the trace 4— >5 overlaps with 
2— >3. The strong maximum in rf diamagnetic response 
and the strong hysteresis at T p are only present during 
the first temperature sweep after ramping the magnetic 
field, after which it only manifests as a kink in x(T"). 
Therefore, we can conclude that the state resulting after 
ZFC to base temperature and increasing magnetic field 
has a much weaker diamagnetic susceptibility then the 
subsequent warming and cooling the sample in applied 
magnetic field (FC). 

Although, the rf Campbell penetration depth cannot 
be used to rule out the existence of the Bragg glass in 
MgB2, we propose as a possible explanation for the hys- 
teresis associated with the peak effect a model based on 
the idea suggested in Ref. [T5] for an individual pinning 
potential, and extend it to the collective pinning regime. 
Within the collective pinning theory 17 , following nota- 
tions from Ref. [29 , for a density of pinning centers, n, 
the free energy density of the vortex bundle (FLL) within 
the correlation volume, V c , can be written as the sum of 
the pinning and the elastic energy: 
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where Cqq is the shear modulus and C44 is the tilt mod- 
ulus of the FLL, f is the pinning force, £-the correlation 
length, and R c and L c describe average distances over 
which the FLL is distorted due to shear and tilt, respec- 
tively. As detailed knowledge of pinning in MgB2 is not 
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FIG. 5: History dependence of x{T) in MgB1.904C0.096 (main 
figure) and MgB1.85C0.15 (lower inset) for applied magnetic 
field H=35 kOe, following a sequence explained in the text 
and sketched in the upper inset. 



yet available, the inset of Fig. f6^a) shows a qualitative 
plot for arbitrary values of the constants in Eq. [3] Dur- 
ing the ramping of the magnetic field after ZFC the sam- 
ple, within the Bean model, there will be a macroscopic 
distribution B(x) of magnetic field inside the sample, ac- 
companied by a current distribution ^i j = dB(x)/dx 
(inset of Fig. |6jb)). These currents exert a Lorentz force 
(/l = J'Po) on vortices causing additional bending of 
the vortex bundle from its equilibrium (minimum en- 
ergy) configuration. The displacement takes place until 
the Lorentz force is compensated by the pinning forces 
(dV/dr = j(j>o)- The situation is sketched in Fig. ffl a) 
and (b). At the peak effect, whether because the pinning 
becomes strong enough or because the elastic modulus 
softens, the vortices relax to equilibrium. The Campbell 
penetration depth will then map the curvature of the 
pinning potential during the relaxation process. As seen 
in Fig. ISFa) , the curvature can change sign, so that the 
absolute value shows a small kink, very similar to that 
observed in the penetration depth of the undoped or Li- 
doped MgB2. If in Eq. [3] the number of pinning sites, 
n is increased (Fig. plb)) while all the other parameters 
are kept the same, the kink becomes more pronounced, 
and the behavior is very similar to that observed in our 
C-doped MgB2 data. Therefore, the model suggests that 
substitution with C at B sites, adds disorder in B planes, 
increasing the number of pinning centers. On the other 
hand, the fact that there is no clear change in the shape 
of the anomaly for Li substitution is consistent with the 
idea that Li does not increase the density of pinning sites. 
Indeed, it slightly lowers the pinning site density as was 
inferred from Fig. HI 

When the sample is further FC, the field distribution 
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FIG. 6: The pinning potential V(r) and the second deriva- 
tive \d 2 V(r)/dr 2 \ (Labusch parameter) of a FLL for arbitrary 
fixed values of Cee, C44, £, f in Eq.|3] and two different values 
for n: ni (a) and 3ni(b). In both graphs it is sketched the 
additional displacement due to the Lorentz force produced by 
screening supercurrents j. Inset (a): The total free energy 
of a vortex bundle as a function of distortion length and the 
contribution from elastic (first two terms of Eq. [3| and pin- 
ning forces (last term of Eq. pi). Inset (b): Field and current 
distribution within the Bean model when the field is ramped 
after ZFC (continuous lines) and after FC (dashed lines). 



inside is uniform (dashed lines of the inset in Fig. p[b)), 
the supercurrent j = and a stable FLL is established. 
This explains both the absence of a clear dip in rf sus- 
ceptibility and the absence of hysteresis for subsequent 
warming and cooling under applied magnetic field. Then, 
the Campbell penetration depth features only a kink at 
the peak effect, consistent with a transition from a weak 
(strong) to strong (weak) pinning regime as the sample 
is warmed (cooled) in field. 

It must be noted that our findings are similar to those 
of Refs. [12)13] . but in previous reports there is still a 
visible minimum in ac susceptibility after FC the sample. 
However, both in Ref. [T^] and Ref. [13], the ac magnetic 
field is at least two orders of magnitude larger than in 
our experimental setup. As explained in the experimental 
section, our technique probes the curvature of the pinning 
potential much more locally (over a few A based on an 



estimate similar to that in Ref. [16 ) whereas in other 
techniques the amplitude of the vortex oscillation is much 
larger. 



IV. CONCLUSIONS 

In conclusion, we have studied the small-amplitude rf 
susceptibility of pure, C, Li and C-Li doped MgB2 single 
crystals. At the peak effect, the rf susceptibility x(T) 
shows a kink, similar to that observed in local ac sus- 
ceptibility measurements. Doping with C was found to 
have a strong influence on the peak effect, consistent 
with increasing the density of pinning centers, whereas 
Li substitution seems to slightly reduce the pinning in 
MgB2- This could have practical implications in find- 
ing the proper way to achieve high critical currents and 
critical magnetic fields in MgLV We propose that the 
history dependence associated with the peak effect is not 



necessarily a signature of a transition to a vortex glass 
phase. Rather it may be the result of additional dis- 
tortion of the FLL caused by screening currents induced 
during a magnetic field ramp. Such a model may have 
broader implications, possibly explaining the history ef- 
fects observed in ac susceptibility measurements on other 
superconductors . 
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